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The interaction of three jets is studied both theoretically and 
experimentally o The angle of deflection and the velocity and turbulence 
distributions were determined through the use of a hot-wire anemometer* 
The free -streamline theory has been applied to the determination of the 
jet deflection and the results so obtained were found to be in good 
agreement with those obtained experimentally o The study is a convincing 
evidence of the fact that a proper combination of inviscid flow analysis 
together with experimental results can provide the necessary guidance 
for the design of proportional amplifiers* 



LI3R 1RY 

N.W/-L T ’3.-rGP-f'TJA'!’E SCHOOL 
IIOMCLHE':, C." -?0 



Library 
Naval Post 
Monterey, 



TABLE OF CONTENTS 



I. INTRODUCTION ---------------------13 

II. THEORETICAL CONSIDERATIONS -------------- 16 

III. EXPERIMENTAL EQUIPMENT AND PROCEDURE -------.-23 

IV. EXPERIMENTAL PROCEDURE ----------------26 

V. DISCUSSION OF RESULTS -----------------27 

VI. CONCLUSIONS ----------------------30 

APPENDIX I EVALUATION OF INTEGRAL (20) ----------- 3 ! 

APPENDIX II EVALUATION OF INTEGRAL (30) -----------37 

APPENDIX III EVALUATION OF INTEGRAL (40) -----------40 

APPENDIX IV EVALUATION OF INTEGRAL (49) -----------43 

APPENDIX V EVALUATION OF DISCONTINUITY OF INTEGRAL (61) - - -45 

APPENDIX VI DERIVATION OF TURBULENCE INTENSITY EQUATION - - - 49 

COMPUTER PROGRAM A ITERATION PROCESS ----------- 117 

COMPUTER PROGRAM B DATA REDUCTION -------------124 

REFERENCES -------------------------- 126 

INITIAL DISTRIBUTION LIST ------- -------- I2 7 

FORM DD 1473 ------------------------- 12 9 



-rf? 

Oc-iilo 



3 



LIST OF ILLUSTRATIONS 



1 . 


Transformations planes 




52 


2 . 


Jet assembly 






53 


3. 


Jet assembly 






54 


4. 


Calibration equipment 




55 


5. 


General ! 


setup of experiment 




56 


6 . 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, T = 0 


deg, r/a * 


0 . 8 ) 


57 


7. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, £ - 0 


deg, r/a - 


1 . 2 ) 


58 


8 . 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, £ - 0 


deg, r/a - 


2 . 0 ) 


59 


9. 


Velocity 


profile and 


turbulence 


intensity 






(Re * 


8440, S'- 0 


deg, r/a = 


2 . 8 ) 


60 


10 . 


Velocity 


profile and 


turbulence 


intensity 






(Re * 


8440, S' = 0 


deg, r/a - 


4.0) 


61 


11 . 


Velocity 


profile and 


turbulence 


intensity 






(Re s 


8440, S~= 0 


deg, r/a - 


8 . 0 ) 


62 


12 . 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, £ = 0 


deg, r/a - 


16.0) 


63 


13. 


Velocity 


profile and 


turbulence 


intensity 






(Re - 


8440, £ = 6 


deg, r/a ■ 


0 . 8 ) 


64 


14. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440 5 5 = 6 


deg, r/a s 


1 . 2 ) 


65 


15. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, S'- 6 


deg, r/a - 


2 . 0 ) 


66 


16 . 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, 6 


deg, r/a = 


4.0) 


67 


17. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, T = 6 


deg, r/a - 


8 . 0 ) 


68 


18. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440 , q = 0 deg, r/a s 0.8) 


69 



5 



19. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, <3^ 0 


deg, r/a =» 


1.2) 


70 


20. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, <r- 0 


deg, r/a = 


1.6) 


71 


21. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, 0 


deg, r/a = 


2.0) 


72 


22. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, r- 0 


deg, r/a - 


2.8) 


73 


23. 


Velocity 


profile and 


turbulence 


intensity 






(Re - 


8440, %= 0 


deg, r/a - 


4.0) 


74 


24. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, 6 


deg, r/a = 


0.8) 


75 


25. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, T= 6 


deg, r/a = 


1.2) 


76 


26. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, f= 6 


deg, r/a - 


1.6) 


77 


27. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, £=6 


deg, r/a - 


2.0) 


78 


28. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, tT= 6 


deg, r/a = 


2.8) 


79 


29. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, r- 6 


deg, r/a - 


4.0) 


80 


30. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, 12 


deg, r/a - 


0.8) 


81 


31. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, S~= 12 


! deg, r/a s 


» 1.2) 


82 


32. 


Velocity 


profile and 


turbulence 


intensity 






(Re =- 


8440, p* 12 


deg, r/a = 


* 1.6) 


83 


33. 


Velocity 


profile and 


turbulence 


intensity 






(Re - 


8440, S'- 12 


: deg, r/a = 2.0) 


84 


34. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, 12 


: deg, r/a = 


» 2.8) 


85 


35. 


Velocity 


profile and 


turbulence 


intensity 






(Re = 


8440, S' =*12 


deg, r/a = 


• 4.0) 


86 


36. 


Velocity 


profile and 


turbulence 


intensity 






(Re « 


8440, o= 12 deg, r/a - 


' 8.0) 


87 



6 



37. 


Velocity 
(Re = 


profile and 
6330, £ = 0 


turbulence intensity 
deg, r/a * 0.8) 


88 


38. 


Velocity 
(Re = 


profile and 
6330, <T = 0 


turbulence intensity 
deg, r/a = !<>2) 


89 


39. 


Velocity 
(Re = 


profile and 
6330, r- 0 


turbulence intensity 
deg, r/a - 1.6) 


90 


40. 


Velocity 
(Re = 


profile and 
6330, £ = 0 


turbulence intensity 
deg, r/a = 2.0) 


91 


41. 


Velocity 
(Re - 


profile and 
6330, f = 0 


turbulence intensity 
deg, r/a = 2.8) 


92 


42. 


Velocity 
(Re = 


profile and 
6330, S = 0 


turbulence intensity 
deg, r/a = 4.0) 


93 


43. 


Velocity 
(Re = 


profile and 
6330, <T = 6 


turbulence intensity 
deg, r/a = 0.8) 


94 


44. 


Velocity 
(Re =-' 


profile and 
6330, <?* = 6 


turbulence intensity 
deg, r/a - 1.2) 


95 


45. 


Velocity 
(Re = 


profile and 
6330, r = 6 


turbulence intensity 
deg, r/a - 1.6) 


96 


46. 


Velocity 
(Re = 


profile and 
6330, f = 6 


turbulence intensity 
deg, r/a =2.0) 


97 


47. 


Velocity 
(Re = 


profile and 
6330, £■= 6 


turbulence intensity 
deg, r/a = 2.8) 


98 


48. 


Velocity 
(Re - 


profile and 
6330, 6* = 6 


turbulence intensity 
deg , r/a = 4.0) 


99 


49. 


Velocity 

(Re-= 


profile and turbulence intensity 
6330, 12 deg, r/a - 0.8) 


100 


50. 


Velocity 
(Re = 


profile and turbulence intensity 
6330, c5"= 12 deg, r/a - 1.2) 


101 


51. 


Velocity 
(Re - 


profile and turbulence intensity 
6330, C - 12 deg, r/a = 1.6) 


102 


52. 


Velocity 
(Re = 


profile and turbulence intensity 
6330, & = 12 deg, r/a = 2.0) 


103 


53. 


Velocity 
(Re = 


profile and turbulence intensity 
6330, 12 deg, r/a = 2.8) 


104 


54. 


Velocity 
(Re = 


profile and turbulence intensity 
6330, <T= 12 deg, r/a - 4.0) 


105 



7 



55. Velocity profiles 

(Re = 6330, 7= 0 deg, (VC+VG)/2VA =0.4) 106 

56. Velocity profiles 

(Re = 6330, 6 deg, (VC+VG)/2VA = 0.4) 107 

57. Velocity profiles 

(Re = 6330, $ = 12 deg, (VC+VG)/2VA = 0.4) 108 

58. Turbulence intensity 

(Re = 6330, f= 0 deg, (VC+VG)/2VA = 0.4) 109 

59. Turbulence intensity 

(Re = 6330, S' = 6 deg, (VC+VG)/2VA = 0.4) no 

60. Velocity profiles at different stations downstream m 

61. Deflection angles (VC+VG)/2VA = 1.0, 0.7, 0.4 112 

62. Deflection angles (VC+VG)/2VA = 0.9, 0.6, 0.3 lb 

63. Deflection angles (VC+VG)/2VA = 0.8, 0.5, 0.2 114 

64. Experimental and theoretical deflection angles 115 

65. Deflection angle for different pressures average 116 



8 



LIST OF SYMBOLS 



b 

CF (P) 
d 

3 

i 

PA 

PC,PG 

q 

st 

r 

Re 

s 

t 

U 

U 




VA 

VC 

VG 

VRO 

VRI 

Vj 



Power nozzle width 
Control nozzles width 

Pressure correction factor for rotameter 
Resultant jet's width at vena contracta 
Imaginary part of the expression 



Power jet, stagnation pressure 
Control jets, stagnation pressure 
Magnitude of local velocity 
Real part of the expression 
Distance from pivoting point 
Reynolds number 
Setback of control nozzles 
Coordinate in t~plane 
Instantaneous velocity 
Local mean velocity 

RoMoSo value of the fluctuating component of velocity 

Power jet velocity 

Left control jet velocity 

Right control jet velocity 

Reference velocity VRO s 0«3 VA 

Reference velocity VRI = 0.4 VA 

Resultant jet velocity 




9 



s 

a 



Complement of the deflection angle 
Deflection angle of resultant jet 

Coordinate in the -plane 



10 



ACKNOWLEDGEMENT 



The work described herein was made possible through the sponsorship 
of the U. S. Army Research Office in Durham, North Carolina. I wish 
to express my appreciation to Dr. T. Sarpkaya for his inspired guidance 
and advice during the course of this investigation. A special note of 
appreciation is also extended to my wife for her patient encouragement 
during this work. 



11 



I 




I . INTRODUCTION 



The interaction of compressible or incompressible jets has been 
subject to study for at least the past fifty years for the purpose of 
predicting the velocity and turbulence distributions in the combined 
jets and in general for the purpose of predicting the performance char- 
acteristics of all sorts of machinery where such jet interaction plays 
a significant role. Most of the theoretical analyses have been restricted 
to the study of the interaction of jets of an inviscid fluid , and the 
complex effects of viscosity and the effects of the bounding walls have 
been largely ignored. 

The advance of the fluidics technology and the discovery, among 
others, of proportional amplifiers have renewed a strong interest in 
both the theoretical and experimental investigation of jet interaction. 
Proportional amplifiers were designed and operated either through the 
use of the "pressure control" concept or through the use of "momentum 
control" concept. In either case, both the single-sided and double- 
sided devices were considered. The pressure controlled device is 
operated on the basis of the differential pressure between the two 
cavities placed on the two sides of the resultant jet. The vents are 
not open to the atmosphere. The momentum controlled proportional ampli- 
fier operates, as its name implies, on the momentum balance of the com- 
bining jets. Thus, it is apparent that, whereas the pressure-controlled 
amplifier depends primarily on pressure to deflect the jet, the momentum- 
controlled amplifier depends primarily on momentum. For certain con- 
figurations where the venting of the side cavities is not complete or 
where the entrainment from the two sides of the jet is not completely 
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attainable, then it is necessary to take both pressure and momentum ef- 
fects into account in attempting to design amplifiers optimized for 
particular conditions. It should, in passing, be noted that the use 
of both control mechanisms tends to produce devices whose band-widths 
are relatively small, even though they may have significantly higher 
pressure gains. We shall confine our attention here to a thorough study 
of the momentum-controlled devices or beam-deflection amplifiers, both 
for the purpose of understanding their fluid mechanics and for the pur- 
pose of contributing to a predictable design of proportional amplifiers. 

The literature concerning the fundamental behavior of beam-deflec- 
tion devices is still relatively limited but is growing rapidly with the 
increasing realization of the design limitations of current fluidics 
technology. 

The prediction of the performance of such devices depends on a 
knowledge of the exact mode of interaction of the jets. There are, 
obviously, many possible geometric and kinematic variables for the 
interaction process and it is not economically feasible to experimental- 
ly investigate the characteristics of the mode of operation of each and 
every possible combination. This realization necessarily led, in the 
first decade of fluidics, to the evolution of several flow models with 
varying degrees of approximations and complexities. There are, to be 
sure, essential differences, as well as similarities, among these models, 
and all must be carefully assessed before a design choice is made. 

Dexter [Ref. 1] employed the "’submerged jet”” model, assuming that 
viscous mixing with the surrounding fluid has reached the stage where 
the velocity profile of the deflected jet is of Gaussian form. Moynihan 
and Reilly [Ref. 2] employed the "’self preservation" 9 model in which each 
jet is assumed to have retained its identity during interaction and is 
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not subjected to any viscous effects as it proceeds downstream 0 Douglas 
and Neve [Ref* 3] extended the works in the references just cited by 
applying the two models to the appropriate regions of the resultant 
jet * 

The present work presents the theoretical results obtained through 
the use of the free -streamline theory for symmetric control-jet con- 
figurations and the experimental results obtained with a hot-wire 
anemometer o 
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II. THEORETICAL CONSIDERATIONS 



As in all applications of the f ree-streamline theory, the pressure 
and magnitude of velocity along the free-stream lines are assumed con- 
stant. It is further assumed that both the pressure and velocity are 
continuous along the separation lines between the power jet and the 
control jets. 

Experiments show that each control jet retains a large measure of 
its integrity, and the power jet undergoes an acceleration or a "vena 
contracta" effect during the initial face of the interaction. The 
combined jet attains a Gaussian profile only after a distance of two 
to six nozzle widths from the pivoting point, depending on the flow 
ratio and the geometry of the system. 

It is apparent from the foregoing that the flow in the region in 
which the momentum-exchange and the turbulent mixing take place does 
not, strictly speaking, follow the assumptions made in the free-stream- 
line theory. The complicated interaction between the jets, and the 
pressure of turbulent mixing, make the analytical predictions of 
three-dimensional flow field extremely difficult. To obtain analytical 
solutions, a simplified incompressible flow model, which retains some 
of the major features of the actual flow field, must be employed. Al- 
though the probable error introduced through the use of such a model is 
very difficult to estimate, previous analyses place a reasonable degree 
of confidence in the applicability of the model to certain features of 
the jet interaction phenomenon. 

A. TRANSFORMATION PLANES 

A schematic drawing of the flow configuration, together with the 
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resultant jet, is shown in Fig* la* This configuration may be trans- 
formed to an-TL-plane through the use of the Planck 0 s transformation 
with the usual assumptions of the f ree-streamline theory . Thus, writing 



XT Xi t is 

and carrying out the necessary manipulations, one obtains Fig. lb. 



( 1 ) 



It is a well-known fact that any polygon, such as the one shown in 
Fig. lb, may be transformed to the upper or the lower half of a real 
plane through the use of Schwartz-Christof fe 1 transformation given by 






^ d t- 



(t - ( t Ct-Cr)^ 



-I si (2) 



where M and N are constants to be determined and Z^ and are given 

by 7 

C, = o 



Then one has 



^ = l 

t/x - 

2 . 



^ = k 



= It % =■ x, 

2 . 



Cr“- W. 

<r= X, 



XI = 



tdfc - 






+ M 



(3) 



The evaluation of this integral is shown in detail in Appendix I. 
The result is that 



.a* jjbg -for t*-k 






(4) 
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Table I gives the corresponding values between the points in the -plane 
and the t -plane shown in Fig. lc. 

The relationships between the parameters shown in the t-plane and 
various velocity ratios are obtained as follows: 

Through the use of the conditions at point A 



_ -la VA 



v i 



-f Ky 5 






one has 



n _ Vft _ * 

( Yi i -v / » - k 2. 

Through the use of the conditions at point C 

-tn^ + lo = v( JEEEjUEE) 

Vj V \l c*-ic» / 



one has 



Of* VC = 



7 C z - k 1 ' 



v j 7T- + V »-c> 



Through the use of the conditions at point G 



— in - 'U ( a 7 

one has Xi \ - *J^ x -k x ; 



* * T 






7 i-k 1 ' a- *h- ( % x 

Through the use of the conditions at point E 



one has 



6 = W ( V>-V 
1 v ; 



sm ft « 



r-\l^ 



(7) 



( 8 ) 



(9) 



( 10 ) 
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Finally, through the use of the conditions on the free streamline, we 
obtain 

-a.- lb . l ±h-*' 4 /r? 



or 



e i& - ± ,nm 

^ e-~ 



4 






fc'-k*- 



t*-k 2 



xj. 



£ - OOS. © 4- i surv © 



(11) 



It should be noted that the conservation of mass yields 



d ~ 4 . b + i^ol 



The t-plane is comprised of three sources and one sink. Thus, the 
complex potential for the flow in the upper half of the t-plane may be 
written in terms of the known potentials of sources and sinks as: 

vV ss ^l(t b\ t _ V C ^ + ^ Q 2 ) 

where W - ^ and (f) is the complex potential and ^ 

the stream function. 

Taking the derivative of w with respect to t, making use of Eqs. 

(7), (8), and (9), and the conservation of mass relationship, one has 



dw _ V; _ d + ^ a. + ^ b + # b 

cJt TL L t-kf t t-VC. t-% J 



(13) 



The distances between various points in the z^plane, the contraction 
coefficient, and the angle of deflection may be evaluated as follows: 



DP 



= a 



J dz j 



(14) 
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and noting that DF = a + 2s, one has 
OL+ ZS 



■ «• [ f * &) dt J 



dw /A 

Replacing — from Eq. (13) and SL from Eqs, (4) and (6), and taking 
the Cauchy principal value of the integral to account for the discon- 
tinuity at t = -f 9 we finally have* 

,-K 

P ' ' 



CL + 2.5 - d 6in5 + 

£-*0 



It-fc *- ( ^ _dC_ 

L 1 i V t t+c 



iwv 



£-*•0 
Q -*oo 




It 

t 



^ ^ b _ d \ dt 

t+c b + jJf&T 



(15) 



lv*v\ 

Q-* oo 
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Juk 1 




Jr lb * _ 


_d_ 


I « 


X - 


v t 
Q 


t+C 




' J 



The details of the integration are shown in Appendix II. The results 
may be expressed as; 



QL. 4 2LS 
lo 



d + 2 h-^ \u 'hz£ 

r x L * b r\ r^- + ^ 



( 16 ) 



_ JL_ U 77 ‘ U\ — — 

*I\-V Z -V k*- ^ N k" \/ »- fc* * \i ■$ - , * t - 

Equation (16) establishes the relationship between the various jet 
velocities and the amplifier setbacko Another expression is needed to 
establish the symmetry of the geometrical conf iguration e This may be 
done by performing an integration from D to F along the free-streamline 
and writing the imaginary part equal to zero so as to ensure that the 
two sides of the control jets are along the same line and that there is 
no unsymmetry. Thus, writing 
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(17) 



cJ 



J d aj= o 



or 






^ o 



and replacing — from Eq. (13) and from Eqs* (4) and (6 ) 9 one has 






-1 



{kk M ^ 1 4. *k _L_ + Ik _d_ - d _J_ ^ dt = o 
x t x, -t-vc. x / 



(18) 



Details of the integration are shown in Appendix III. The analysis 
yields: 



X O + O 




iiiliW +L^]y - 
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V + 4 W \ d_ 
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(k4»')(cvi > ) ) 


k.4 \ I 
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2(k + {) 


2^\ 


(M »)({-m) 


k+i / 



- 0 



(19) 
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where K, C^) is the complete elliptic integral of the first kind and 

K j is the complete elliptic integral of the third kind* 

2 

Since the parameter^ of the elliptic integral of the third kind 

2 

shown in Eq. (19) was not always in the range of 0 1» a sub- 

routine program was written to reduce all the complete elliptic inte- 
grals of the third kind into the sum of complete and incomplete elliptic 
integrals of the first and second kind for which there are tables and 
subroutines available at the computer facility of the Naval Postgraduate 
School. 

The reductions were made through the use of the transformations 
given in Ref* [4], and through the use of Heuman's Lambda and Jacobi's 
Zeta functions, as given in Ref* [5]. The aforementioned reductions are 
shown in Appendix IV. The evaluations of the resulting expressions for 
the discontinuity at t = -f is shown in Appendix V. The results obtained 
for - \ and b x | through the use of the expressions given above 

a. cC 

will be presented together with the experimental data in connection with 
the discussion of results* 
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The experimental apparatus consisted of a jet assembly*, hot-wire 
aneomometer system*, and a velocity calibrator 

A. JET ASSEMBLY 

As shown in Fig, 2 and 3*, the entire assembly consisted of three 
nozzles, one for the power jet and two for the control jets. Each noz- 
zle was 1/4 inch wide and 2 inches high. This resulted in an aspect 
ratio of eight. The jet assembly was sandwiched between two 1/4 inch 
plexiglass plates. 

The width of the control nozzle was adjustable from 1/4 inch to 
1/2 inch. For all the tests reported here, the control and power jet 
nozzles, as well as the setback, were maintained at widths of 1/4 inch. 

The smooth streamline contour of the nozzles gave an essentially 
uniform exit velocity profile. The special configuration adopted for 
the control jet walls prevented the jet from attaching to the side walls 
through the Coanda effect. In addition, two large holes, as shown in 
Fig, 2, were provided in the upper and lower cover-plates, in order to 
minimize the restriction to entrainment to the resultant jet. 

The supply lines were 1-1/4 inch diameter Tygon tubing, connected 
to the three inlets. The flow was straightened via a one inch long 
honey-comb section before undergoing a contraction of ten to one. 

Supply pressures were measured by means of inclined manometers 
through the use of wall pressure tabs located in the straight part of 
each nozzle. 

The top and bottom plates had two slots cut out to enable the hot- 
wire anemometer mechanism to traverse along lines parallel to the power 
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nozzle axis* A graph paper * graduated in degrees * was placed under the 
bottom plate for the purpose of quickly measuring the jet deflection 
angle by means of a small light string attached to a needle* A more 
precise measurement of the jet deflection angle was* of course* made 
by means of the hot-wire anemometer* 

Air was supplied to the system from a 50hp compressor at approxi- 
mately 160 psi. Air first passed through a micro-filter to remove oil 
and dust and then through the primary pressure regulator* The three 
jets were independently supplied from this regulator* Each jet supply 
line had on it a series of pressure regulators, thus providing indepen- 
dently regulated flows into three rotameters* The rotameters were 
identical and each had a maximum capacity of 17*5 scfm* at 14*7 psia 
and 70°F. The actual flow rates were obtained from 

Q = (rotameter reading* % of full flow) x (17*5) x CF(P) 
where CF(F) represents the pressure correction factor obtained from a 
chart provided by the manufacturer* 

B* HOT-WIRE ANEMOMETER 

A constant temperature hot-wire anemometer* manufactured by the 
Thermo-Systems, Model 1050* was used in the measurement of velocity and 
turbulence intensities* The wire used (wire #1272-T1*5) had a cold re- 
sistance of 7*30 ohms and an operating resistance of 10*95 ohms* The 
voltage across the bridge was measured with a digital voltmeter accurate 
to 0*01 volts* 

A special probe holder was constructed and is shown in Fig* 4* This 
holder which was fitted into the slots in the jet assembly* enabled the 
probe to move along lines parallel to the power-nozzle axis in the mid- 
plane. A micrometer mechanism on the holder provided transverse move- 
ment of the hot-wire across the jet with an accuracy of 0*001 inch* 
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The hot-wire probe was calibrated with the Thermo-Systems Cali- 
brator, Model 1125, shown in Fig. 4. An inclined manometer was con- 
nected to the reservoir of the calibrator to measure the stagnation 
pressure and the hot-wire was placed 1/16 inch above the 0.150 inch 
orifice of the calibrator. For the details of calibration analysis, 
see Appendix VI. 
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IV o EXPERIMENTAL PROCEDURE 



The experiments were carried out for a combination of three power 
jet velocities and for four control jet bias velocities. The bias 
velocity was defined as (VC+VG) /2VA. For each combination, the velo- 
city and turbulence distributions were measured at various distances from 
the pivoting point, and the results were reduced to actual data through 
the use of the computerized calibration curve * Most of the experiments 
were repeated for the purpose of checking the reproducibility of the 
data. Finally, the deflection angles were determined from the velocity 
profiles by connecting the peaks of each successive profile with the 
pivoting point. The deflection angles obtained in this manner did not 
deviate more than +0-5 degrees from those obtained through the use of 
the floating string and chart paper combination. 

The supply pressure for each jet was measured, as previously cited, 
by means of inclined manometers. The stagnation pressure of the power 
jet was measured by one manometer, and the stagnation, as well as 
differential pressure of the control jets, through the use of two other 
manometers. 
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V. DISCUSSION OF RESULTS 



In the following, first the experimental results obtained with the 
hot-wire anemometer will be discussed and then a comparison will be 
made between the theoretically and experimentally obtained deflection 
angle. 

Figures 6 through 54 represent the velocity and turbulence intensity 
distributions as a function of x/a and U/VRO in terms of normalized 
control jet velocities, deflection angle, and the distance from the 
pivoting point. It is apparent from all of these figures that the con- 
trol jets retain their identity up to a distance of 2a from the pivoting 
point. This region will be called the mixing region. For distances 
larger than approximately 2a, the velocity profiles acquire a nearly 
Guassian distribution. The turbulence intensity distribution exhibits 
essentially four peaks: two in the outer edges of the resultant jet 
because of the mixing of the jet with the ambient flow and two peaks in 
the regions where the control jets mix with the power jet. As the dis- 
tances increase from the pivoting point, the turbulence intensity in the 
inner regions decreases and eventually reduces to a single minimum at the 
axis of the jet. It is therefore apparent that the application of the 
free-streamline theory to the mixing of jets of unequal velocities is 
not entirely valid in the mixing region since there are three independent 
free-streamline velocities. In the fully established region, the velocity 
profile is Gaussian and thus far from being uniform. Consequently, the 
use of the single jet velocity is not fully justified unless special pre- 
cautions are taken to make sure that the gross continuity of the mass 
flow is satisfied by properly selecting the apparent jet width. In so 
doing, the maximum centerline velocity may be used as the free-streamline 
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velocity and the effective jet velocity may be calculated from the equa- 
tion of continuity . The theoretical analysis presented herein makes 

good use of this fact and shows that the results so obtained are in good 
agreement with those obtained experimentally. 

Figures 55 through 57 show a superposition of the velocity profiles 
and the evolution of the mixing from three independent jets to a single 
profile, Similar graphs are presented in FigSo 58 and 59 for the 
turbulence intensity.. Figure 60 is a three-dimensional plot of the 
evolution of the velocity profiles as a function of the distance from 
the pivoting point., 

The deflection angle is shown in Figs, 61 through 63 as a function 
of the sum and difference of the normalized control jets for all 
Reynolds numbers. It is immediately apparent that the deflection angle 
is a linear function of the normalized differential control jets velocity 
and that the deflection angle does not vary appreciably with the normal- 
ized average of the control jets. These results suggest that the 
characteristics of beam*def lection amplifiers be evaluated in terms of 
the parameters cited above rather than those previously used by other 
investigators. Although the experimental results presented herein are 
for equal control and power jet nozzles, there is every reason to be- 
lieve that similar correlations will exist between the jet deflection 
angle and the two control jets parameters, 

A comparison of the experimental results with those obtained 
theoretically for a normalized average control jet velocity of unity 
is shown in Fig, 64, As previously cited, the final jet width d used 
in the calculations was obtained by dividing the total flow rate by the 
maximum jet velocity. It is believed that there is at present no other 
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way of applying the free-streamline\ theory to the deflection of viscous 
jets for in the mixing region there are three different jet velocities 
and three fairly uniform. velocity distributions and, in the fully estab- 
lished region, a single maximum with a Gaussian distribution. Thus, the 
only logical assumption to be made is that the free-streamline velocity 
used in the analysis is equal to the maximum of the: fully established pro- 
file and that the effective jet width is determined by the equation of 
continuity. In reality, it appears that the maximum centerline velocity 
is somewhat higher than that predicted by the free-streamline theory and 
that the effective jet width is somewhat smaller than that used in the 
analysis. Be that as it may, no second-order corrections were applied 
to the analysis in order to preserve its generality. The extension of 
the analysis to other velocity ratios will be separately undertaken in 
the future. 

The stagnation pressures at supply reservoirs of the three jets 
were determined as previously explained and are plotted in Fig. 65. 

Mean values rather than the actual data points are shown in this graph 
because the data points fell fairly consistently on straight lines. 

It is immediately apparent that the normalized differential control port 
pressures are linear functions of the deflection angle and depend weakly 
on the mean pressure of the control ports. A simple analysis of the 
jet deflection through the use of the momentum balance equation shows 
that the deflection angle so determined is considerably higher than that 
obtained experimentally or theoretically even when the pressure differ- 
ences between various nozzles are taken into consideration. 
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VI, CONCLUSIONS 



The foregoing theoretical and experimental study warrants the 
following conclusions: 

(a) The jets preserve their identity in a region of approxi® 
mately 2a from the pivoting point; 

(b) The velocity profiles become Gaussian beyond the mixing 
region ; 

(c) The free streamline theory predicts the deflection angles 
fairly accurately provided that an effective jet width 
(calculated in terms of the mass flow and the maximum 
centerline velocity) is incorporated into the analysis; 

(d) The previous models cited in the introduction do not 
predict, or are not capable of predicting, the jet de - 
flection angle to an accuracy demanded by the beam de- 
flection amplifiers; and finally, 

(e) The normalized differential control reservoir pressures 
vary linearly with the deflection angle. It is thus 
believed that an analysis such as the one presented here, 
together with one experimentally determined parameter, is 
fully capable of predicting the performance of proportional 
amplifiers. The knowledge of the jet deflection angle and 
the Gaussian shape of the velocity profile provide through 
a simple integration the determination of the amounts of 
flow at a given amplifier receiver. 
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APPENDIX I 



EVALUATION OF THE INTEGRAL 
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The constants and may be determined, for the different regions of 
the t plane, as follows: 
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Finally, 
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APPENDIX II 



\ 



EVALUATION OF THE INTEGRAL 
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Equations (31), (32) and (33) have no discontinuity at t = - f 
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For Eq. (31), one has 
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Replacing Eqs. (35), (36), (37) and (38) in Eq. (15) one has 
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APPENDIX III 



EVALUATION OF THE INTEGRAL 
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Note that the first integral in Eq„ (42) is an odd function of t. Thus 
its integration from -1 to +1 reduces to zero. The same result is 
obtained for the third integral when p « 0. 
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In taking the limit of Eqs. (43) and (44), it must be noted that the in- 
tegral from -1 to +1 is to be performed by going from -1 to - 90 , to + 50 
to +1. Replacing the value of t in Eq. (46), it is seen that ij> goes from 
^ to zero when t goes from -1 to + 1 , in the sense described above. 
Replacing Eqs. (43) and (44) in Eq. (42) and noting that 'j' goes from 
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where is the complete elliptic integral of the first kind and 

h is as defined by Eq. (45). 
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Replacing Eq* (47) in Bq. (40) and replacing p by its appropriate ,value g 
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The singularity at t = -f is evaluated separately in Appendix V. 
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APPENDIX IV 



EVALUATION OF COMPLETE ELLIPTIC INTEGRAL OF THE THIRD KIND 



In the general case of Complete Elliptic Integrals of third kind 
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where 
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E ^ ) is elliptic integral of the second kind (57) 
F is the elliptic integral of the first kind (58) 
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Subroutine CEL3 was written to calculate the complete elliptic integrals 
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of the third kind for all the cases of cited above. 
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APPENDIX V 



EVALUATION OF THE DISCONTINUITY OF THE INTEGRAL 
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( 66 ) 
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Simplifying, Eq. (69) reduces to 
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Using L'Hopital/s rule 
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APPENDIX VI 



DERIVATION OF TURBULENCE - INTENSITY EQUATION 



On a hot-wire anemometer the heat transfer equilibrium equation is 
given by 
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Rw - operating resistance 
^ = cold resistance 
A And B - constants 
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in other words, turbulence intensity is equal to 
times RMS reading* 



4 E 
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TABLE I 



Point in the 
z-plane 


Corresponding Value 
in the 
t -plane 


Corresponding Value 
in the 
-plane 


A 


0 


- U w til 

Vi 2. 


B 


-k 


i c 
l Vz 


C 


-c 


- U +ic 


D 


-1 


c + i o 


E 


f 


0 4 i fv 


F 


1 


0 + (- TCy 


G 


g 


- W -y t H 

Yi 

, j* 


H 


k 


+ oc -t i , tt v 
l 'X 

* 
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Figure lb -plane 
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Figure lc t-plane 



Figure 1 Transformations Planes 
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Figure 2 Jet Assembly 
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Figure 3 



Jet Assembly 
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Figure 4 Calibration Equipment 




Figure 5 General Setup of Experiment 
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FIGURE 6 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 7 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 8 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 9 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 10 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 11 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 12 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 13 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 14 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 15 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 16 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 17 VELOCITY PROFILE AND TURBULENCE INTENSITY 



68 




FIGURE 18 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 19 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 20 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 21 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 23 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 24 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 25 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 27 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 28 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 29 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 30 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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s FIGURE 31 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 32 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 33 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 34 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 35 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 36 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 37 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 38 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 39 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 40 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 41 VELOCITY PROFILE AND TURBULENCE INTENSITY 



92 




93 




FIGURE 43 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 44 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 45 VELOCITY PROFILE AND TURBULENCE INTENSITY 



96 





FIGURE 46 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 47 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 48 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 49 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 50 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 51 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 52 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 53 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 54 VELOCITY PROFILE AND TURBULENCE INTENSITY 
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FIGURE 56 VELOCITY PROFILES 
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FIGURE 59 TURBULENCE INTENSITY 
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FIGURE 60 VELOCITY PROFILES AT DIFFERENT STATIONS DOWNSTREAM 
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FIGURE 61 DEFLECTION ANGLES (VC+VG)/2VA 
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FIGURE 62 DEFLECTION ANGLES (VC+VG)/2VA =» 0.9, 0.6, 0.3 
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FIGURE 63 DEFLECTION ANGLES (VC+VG)/2VA = 0.8, 0.5, 0.2 
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FIGURE 64 EXPERIMENTAL AND THEORETICAL DEFLECTION ANGLES 







FIGURE 65 DEFLECTION ANGLE FOR DIFFERENT PRESSURES AVERAGE 
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